Interaction of carbon clusters with Ni(lOO) : 
Application to the nucleation of carbon nanotubes 
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In order to understand the first stages of the nucleation of carbon nanotubes in catalytic pro- 
cesses, we present a tight-binding Monte Carlo study of the stability and cohesive mechanisms of 
different carbon structures deposited on nickel (100) surfaces. Depending on the geometry, we 
obtain contrasted results. On the one hand, the analysis of the local energy distributions of flat 
carbon sheets, demonstrate that dangling bonds remain unsaturated in spite of the presence of the 
metallic catalyst. Their adhesion results from the energy gain of the surface Ni atoms located below 
the carbon nanostructure. On the other hand, carbon caps are stabilized by the presence of carbon 
atoms occupying the hollow sites of the fee nickel structure suggesting the saturation of the dangling 
bonds. 
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I. INTRODUCTION 



During the past few years, there has been a tremendous 
increase in research activities devoted to single- wall car- 
bon nanotubes (SWNTs) mainly because of their extraor- 
dinary mechanical, electrical, and chemical properties. 
Future nanoscale applications will however necessitate 
the production of pure SWNTs in significant quantities 
and also the selective formation of nanotubes with well- 
defined chirality or diameter. This obviously requires to 
understand their growth mechanisms. SWNTs are gen- 
erally produced via techniques working at high tempera- 
tures (arc discharge laser ablation or at medium 
temperature (chemical vapor deposition [3]). The pres- 
ence of a transition metal catalyst (such as Fe, Co, Ni, 
. . . ) has been found necessary for the synthesis of SWNTs 
but the exact role played by the catalyst is still unclear. 
The similarities between the samples synthesized from 
different techniques suggest a common growth mecha- 
nism and the so-called vapor-liquid-solid (VLS) model is 
frequently mentioned 0, In this scheme, the growth 
of SWNTs begins with the formation of a metallic parti- 
cle supersaturated with carbon; then carbon segregates 
towards the surface and finally carbonaceous structures 
emerge from the catalyst surface. From the theoreti- 
cal side, several modelizations at an atomic scale have 
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been performed using either phenomenological potentials 
SSBH or first principles methods [1, H, [ll [3 • 

To some extent these computer simulation works con- 
firm the VLS model but with some differences between 
the high temperature and medium temperature routes. 
In chemical vapor deposition methods the size of the 
tubes is generally determined by the size of the parti- 
cle. Actually several modelizations show that carbon 
atoms at the surface self-organize and form graphene-like 
sheets wrapping the particle. Once formed these sheets 
no longer interact with metallic atoms and can detach to 
form nanotubes. This accounts for a so-called tangen- 
tial process. In other situations, principally within the 
high temperature route, bundles of nanotubes grow per- 
pendicularly to the surface and it is still a challenge to 
understand why this can be more advantageous than a 
tangential growth. 

Some arguments based on classical nucleation and 
growth thermodynamic models have already been put 
forward jT3, [ll]. Carbon atoms at the surface of the 
metallic catalyst are assumed to condense in the form of 
graphene flakes. Then the metallic substrate can help 
to saturate the dangling bonds and this is presumably 
favoured by the formation of a cap, the energy cost due 
to the curvature induced by the presence of pentagons, 
being more than compensated by the reduction of the 
number of dangling bonds il8il . This model has been 
substantiated by Fan et al. |14| who performed ah initio 
energy calculations of different arrangements of carbon 
atoms on a Ni(lOO) surface. Among the different carbon 
geometries tested, they found actually that the most sta- 
ble one is a hemispherical cap forming a tube embryo. 
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These calculations are quite accurate but computation- 
ally very demanding and the structures cannot be fully 
relaxed. Furthermore these calculations do not provide 
an easy access to a local analysis of the bonding. Fan 
et al. have thus developed a simplified phcnomcnological 
analysis of the energies to characterize this effect of dan- 
gling bond saturation ^loj . A more precise analysis of the 
chirality fixed by the embryos growing on a (111) surface 
has been made by Reich et al. 16] based on the use of 
the SIESTA ab initio code. They confirm the bond satu- 
ration effect and show that the possible chiral selectivity 
depends more on the carbon-metal bonds than on the 
energy of the different carbon caps. In this calculation 
however, the metallic (nickel) atoms remain unrelaxed, 
a strong approximation, considering the order of magni- 
tude of the NiC interaction, and no local energy analysis 
can either be achieved. 

In this paper, we study the adhesive properties of car- 
bon clusters on a Ni(lOO) surface using a tight-binding 
(TB) total energy model. This model relies on local 
(atomic) energy calculations using the recursion method 
and is coupled with Monte Carlo (MC) simulations in or- 
der to relax the structures. Our objective is to establish if 
the nucleation of SWNTs is favoured by the elimination 
of dangling bonds, as suggested by the above mentioned 
authors and also to discuss the role of the displacements 
of the metallic atoms. We thus present an analysis of the 
local energy distributions in order to determine which 
atoms (carbon or nickel), and to what extent, are sta- 
bilized when various carbon clusters are put in contact 
with a metallic surface. 



II. METHOD 

In our study, the interaction between metal (here 
nickel) and carbon atoms is treated within the semi- 
empirical tight-binding model described in Ref. [20j . 
Only s, p electrons of carbon and d electrons of nickel 
are taken into account. The total energy of an atom i 
is splitted in two parts, a band structure term that de- 
scribes the formation of an energy band when atoms are 
put together and a repulsive term that empirically ac- 
counts for the ionic and electronic repulsions : 

^tot — ^band + ^rep- (1) 

The total energy of the system, Etot , then writes : 

Etot = ^ El^f (2) 

i atoms 

The band energy, El^^^, is given by : 

ELnd^ r\E-4)n,{E)dE, (3) 

J —oo 

where Ef denotes the Fermi level and e° is the atomic 
energy level. We use the recursion method to calculate 



the local electronic density of states ni{E) on each site 
(23| . Only the first four continued fraction coefficients, 
(oi, 6i, a2, &2) corresponding to the first four moments 
of the local density of states are calculated exactly. The 
continued fraction is then expanded to the N*^ using 
constant coefficient equal to 02 and 62 (here N = 40). 
The local density of states ni{E) is then expressed as a 
set of N poles {El , j = 1, TV) with weights (A^ , j = 1, TV) 
by diagonalizing a tridiagonal matrix of size N . Thus, the 
band energy of an atom i writes : 

Eland ^'Y.'^^I (4) 

In the above equation the highest occupied energy level 
jmax depends on each site i. The set of C-C, Ni-Ni 
and Ni-C tight-binding interaction parameters used was 
shown to be transferable to different atomic configura- 
tions by comparison with ab initio calculations. It re- 
produces the expected phase separation tendency for the 
Ni-C alloy, with an enthalpy of mixing of -1-1.0 eV/atom 
for NiC in the (metastable) rocksalt structure fitted to 
ab initio data (ABINIT code). The heat of solution of 
a C atom in a Ni matrix is then found equal to -1-0.4 
eV/atom, in agreement with the experimental value of 
Shelton et al. (0.49 eV/atom [ll]) and the first prin- 
ciples calculations of Siegel et al. (-1-0.2-0.35 eV/atom 
[2^). More details are given elsewhere [2^. This ap- 
proach which correctly describes the chemical bonds of a 
mixed metal-carbon system, is able to deal with systems 
of a few hundreds of atoms over large time scales. As 
compared to the standard diagonalization technique, the 
recursion method is local and then gives direct access to 
the total energy of each atom i (see Eq|3]). This possi- 
bility to easily analyze local energy distributions will be 
very useful in the following. 

This atomic interaction model is then implemented in 
a Monte Carlo code, based on the Metropolis algorithm, 
which allows to relax the structures [l^] ■ Finally different 
tests were performed to check the reliability of our model. 
As an example, a simple relaxation of C atoms cover- 
ing a Ni(lOO) surface induced the "clock" reconstruction 
observed experimentally by Klink et al. [l3|. This re- 
construction has also been studied theoretically [2^, [l^] 
and corresponds to the rotation of the four Ni atoms sur- 
rounding a C atom in an Ni(lOO) surface, clockwise and 
counterclockwise alternatively (see Fig. [1]) 

In a previous study, we presented Grand Canonical 
Monte Carlo simulations of the dissolution and segrega- 
tion of C in the presence of a Ni(lll) surface [21| in 
order to understand the processes involved in the cat- 
alytic growth of carbon nanotubes. At finite temperature 
(1000 and 1500 K), we analyzed the self-organization of 
carbon atoms close to a nickel surface and identified a 
succession of four types of configurations, corresponding 
to an increase of the carbon chemical potential (and con- 
sequently an increase of the number of carbon atoms): 
single C atoms adsorbed on the surface or incorporated 
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FIG. 1: (a) Top view of a surface coverage equal to 0.5 ML. 
(b) Clock reconstruction leading to a (2 x 2)p4(; symmetry. 



in interstitial sites, chains creeping on the surface, de- 
tached sp^ C layers and finally a 3D amorphous C phase. 




FIG. 2: Stability at 500 K of carbon clusters (top and side 
views) deposited on Ni(lOO). (a) flake of 24 atoms containing 
7 hexagons; (b) flake of 54 atoms containing 19 hexagons; (c) 
curvature of the flake of 20 atoms containing a pentagon; (d) 
embryo of (6,6) nanotube. 



III. RESULTS 

Here, we use Monte Carlo simulations in the canonical 
ensemble (with a fixed number of Ni and C atoms) to 
relax the atomic structures and analyze the stability of 
specific carbon structures deposited on a nickel surface. 
We perform MC simulations on a Ni slab presenting a 
(100) surface. This is a four-layer slab of 288 atoms, 
with bond lengths of 2.48 A, separated by a 20 A thick 
vacuum region along the z axis and periodic boundary 
conditions. The slab size is 17.60 x 19.36 x 5.28 A^. 
We consider several configurations of carbon atoms, with 
bond lengths of 1.42 A, deposited on the metallic surface: 
molecules (C5 and Cg rings), a finite linear chain of 12 
atoms parallel to the surface, flat flakes of 24 and 54 
atoms (see Fig. H^-b) and curved structures such as a 



flake of 20 atoms containing a pentagon (see Fig. [2}:) and 
the cap of a (6, 6) tube (see Fig.[21i). 

First, we study the stability of these carbon clusters 
in the presence of the metallic substrate at 500 K. At 
this temperature, the C5 and Ce rings are dissociated to 
present C-C bonds of roug hly 2.2A and C-Ni bonds of 1.9 
A. The twofold coordinated atoms belonging to the finite 
linear chain interact moderately with the metallic surface 
(in the 1.0 eV/atom range), while we observe that end of 
chain, one-fold coordinated, C atom tends to penetrate 
the metallic substrate, to occupy a semi-octahedral site of 
the Ni(lOO) surface. This process has also been observed 
experimentally by Fujii et al [2^. The stability of the 
clusters seems to depend on the number of of C-C bonds, 
with a large enough number of such bonds resulting in 
a stable cluster. For small number of C-C bonds, the 
competition between the formation of a carbide and the 
phase separation results in the formation of a Ni-C alloy. 
For the other structures, the number of C-C bonds is 
large enough to stabilize the carbon clusters despite the 
presence of the metallic surface. The flakes of 24 (Fig. 
[5^) and 54 (Fig. ^Bp) atoms remain in plane and we also 
do not observe any dissociation for the flake of 20 atoms 
containing a pentagon (Fig. [2}d). The presence of this 
defect induces a curvature of the structure in agreement 
with Eulcr's theorem. Last, the interaction between the 
metallic particle and the open tube cap stabilizes it by 
preventing its closure (Fig.[2ji). 

We then calculate the adhesion energies per carbon 
atom, Ea, given by : 



{Nc + N]S^)Ec/N^ ^ NqEc - N]S^E, 



Ni 



Nc 



(5) 



where Njqi and Nc are the number of Ni and C atoms. 
Ec/Ni, and Ec represent the total energies (per 

atom) of the carbon-covered Ni slab at different surface 
coverages, the clean Ni slab and the relevant isolated 
carbon structure, respectively. They are calculated at 
K on fully relaxed configurations. In each case, the 
minimum energy structures are obtained by performing 
simulated annealing. In this process the initial atomic 
configuration, initially at high temperature (1000 K), is 
slowly cooled down so that the system is at any step 
at thermodynamic equilibrium. The major advantage 
of this method is its ability to avoid becoming trapped 
in local minima. Typical runs consist in 10^ external 
Monte Carlo loops, each of them randomly performing 
10'^ atomic displacements trials. 

For all the structures tested, the negative value of the 
adhesion energies, listed in Table [H means that the car- 
bon cluster is stabilized by the surface. For the flakes 
deposited parallel to the surface, the results clearly show 
that it is favorable to incorporate pentagons during the 
initial stages of nucleation. It is therefore tempting to 
think that the structure containing defects can deform 
into a dome to better saturate the dangling bonds at 
their edges. 

To investigate further, we analyze the local energies of 
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TABLE L Adhesion energies per carbon atom of the carbon structures deposited on the surface of Ni(lOO). 
Configurations Flake Flake Flake Linear chain Capped tube 



(Fig. Hi) (Fig, [la) (Fig.H;) (Fig. Hi) 

Nc 24 54 20 12 60 
Ec (eV/at) -6.67 -6.92 -6.55 -5.99 -6.88 
Ec/Ni (eV/at) -4.32 -4.59 -4.30 -4.22 -4.62 
Ea(eV/at C) -0.27 -0.27 -0.58 -LOl -0.29 



carbon and nickel atoms in all the previous cases in order 
to determine which C and Ni atoms atoms gain energy. 
As an example, we plot our results for the flake of 24 
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FIG. 3: Analysis of local energies for a flake of 24 carbon 
atoms in interaction with a Ni(lOO) surface, (a) Histogram of 
C atoms, (b) Histogram of Ni atoms. 

atoms in a histogram form in Fig. [3] After relaxation 
of the isolated carbon flake (away from the metal sur- 
face) , we identify three populations of carbon atoms (full 
lines in Fig.jS^). First, the three-fold coordinated atoms 
denoted by A : 6 atoms in a local configuration equiva- 
lent to an infinite graphene sheet in which each atom is 
in a three-fold coordinated state with sp^ hybrids form- 
ing strong covalent bonds. The energies lie between -7.5 
and -7.4 eV/atom, close to the total energy of a perfect 



graphene sheet calculated with our TB model. Then, the 
6 three-fold coordinated atoms located at the periphery 
of the structure and denoted B: their energies lie between 
-7.3 and -7.2 eV/atom. The difference between the two 
previous populations is simply due to the changes in their 
second neighbor shell. Finally, we have 12 two- fold coor- 
dinated atoms at the periphery of the flake, with energies 
between -5.8 and -5.7 eV/at. These atoms only have two 
neighbors, so that, in principle, one sp^ hybrid is left 
unbonded, forming a dangling bond. With this analy- 
sis, we can estimate the energy of a dangling bond at 
about 1.80 eV: it is the energy difference, appearing on 
the histogram, between a C atom located at the surface 
and a three-fold coordinated A atom. This is close to the 
1.78 eV [29|] and the 2.26 eV ^M] found in finite cluster by 
LDA studies and smaller than the 2.2 to 3.3 eV [s^l found 
in previous LDA works in the case of a tube. Studying 
now the flake in contact with the metallic surface, we can 
see that the three populations are still present, but mod- 
ified. In Fig. |3|the histogram still presents three clearly 
identified regions (dotted lines) equivalent to those de- 
scribed previously, with a larger distribution due to the 
accuracy of the relaxation and to the fact carbon atoms 
that are equivalent in an isolated flake may become no 
longer equivalent when the flake is deposited on a surface. 
The small differences observed however suggest that car- 
bon atoms remains essentially unaffected when deposited 
on the nickel surface, implying that, in this sense, the 
interaction between the metallic surface and the carbon 
structure do not saturate the dangling bonds. This obser- 
vation is somewhat different from the conclusions drawn 
by Fan et at, based on ab initio calculations, and there- 
fore requires further analysis. In agreement with them, 
we nevertheless notice that the energy balance presented 
in Table HI reveals a weak adhesion energy, meaning that 
the flake is more stable on the surface than free standing. 
Since this stabilization docs not come from the carbon 
atoms, these results imply that nickel atoms should be 
stabilized by the presence of carbon. We confirm this 
with the analysis of the local energies of Ni atoms shown 
in Fig. |5Jd and Fig. m In the case of the bare relaxed 
Ni slab, we observe two different populations (full lines 
in Fig. Os). The first corresponds to the bulk atoms 
which have 12 neighbors and energies between -4.4 and 
-4.3 eV/atom. The total energy of a face centered cubic 
bulk Ni is El^^ =-4.44 cV/atom in our TB model. The 
second population corresponds to the 8-fold coordinated 
surface atoms. These are less stable than atoms in bulk 
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FIG. 4: Histogram of Ni atoms close to the flake of graphene. 
Reconstruction of the Ni(lOO) surface when C atoms are 
present. 



and have energies between -4.0 and -3.8 eV/atom. In 
the presence of the carbon structure, both populations 
are still present. It is important to note that 288 atoms 
of Ni have been reported on this histogram. Thus, the 
difference observed on Fig. [3]d which seems weak when 
compared to the analysis of C atoms where only 24 atoms 
appear, corresponds to a larger contribution to the total 
energy. The quantitative analysis of the local energies 
shown in Fig. 01 demonstrates that the nickel atoms, lo- 
cated at the surface and in direct interaction with the 
carbon structure gain energy. On the one hand, the av- 
erage total energy of nickel atoms at the surface and far 
away from carbon atoms is around -3.90 eV/atom. This 
value is similar to the one obtained without carbon atoms 
(full lines in Fig.[3]3); on the other hand, the nickel atoms 
close to the carbon structure are stabilized (Fig.|3|). The 
energy gain is approximatively equal to 0.4 eV/atom and 
is associated with a slight reconstruction of the Ni sur- 
face underlying the flake, as shown in Fig. U) In the case 
of the other flat structure containing 54 carbon atoms, 
the local energy analysis leads to the same conclusions: 
the local energies of carbon atoms are not modified but 
surface nickel atoms below the flake are stabilized with 
roughly the same energy gain. 

The same energy balance is obtained for the flake con- 
taining a pentagon. For the linear chain, the one-fold co- 
ordinated, end of chain, carbon atoms interact strongly 
with the surface and penetrate into the metallic substrate 
to occupy a semi-octahedral site of the (100) Ni structure 
(Fig. [5]) with an energy gain of 0.5 eV. The same effect ap- 
pears with the tube embryo where some carbon atoms are 
located in the fourfold hollow sites. These results show 
clearly that the anchoring of carbon atoms to the catalyst 
is energetically favorable during the nucleation process. 
They gain energy by occupying the semi-octahedral sites 
of the fee nickel (100) surface. 

Finally the adhesion process of carbon sheets on Ni 
(100) is slightly more complex than anticipated. The 
adhesion energy of flat sheets is mainly due to the energy 
gain of the nickel atoms below these sheets. When they 
curve to form caps, the energy gain becomes concentrated 
on the carbon and nickel atoms atoms close to their edge. 



In this case one might argue that dangling bonds are 
saturated, but finally, this is a weak effect which does 
not play very much in favour of curved caps: the energy 
of adhesion of flat and curved sheets are similar and small 
compared to dangling bond energies, for small clusters at 
least. 

At this point it is useful to stress that it can be difficult 
to choose between various decompositions and normal- 
izations of energies of formation or adhesion. Depending 
on the point of view, it can be fruitful to count ener- 
gies per carbon atom of the whole structure (flat sheet) 
or per edge atom (curved sheet). Our model has the 
advantage to provide us with unambiguous local (site) 
energies. Bond energy analyses could be performed as 
well but would be more complex to use in the present 
context. 




top layer 



^^^j bottom layer 




FIG. 5: lUustration of a C atom occupying a hollow site of 
Ni(lOO) surface. 



IV. DISCUSSION 

Our results allow us more generally to understand 
the mechanisms by which carbon atoms adsorb on the 
Ni(lOO) surface. Experimentally, investigations including 
low-energy electron diffracti on jlOl| . surface extended-x- 
ray-absorption fine structure [ll| , and scanning tunneling 
microscopy jl2l | have provided a large amount of informa- 
tion on the interaction of C with Ni(lOO). When carbon 
adsorbs on Ni(lOO) with a coverage smaller than one third 
of a monolayer (ML), it occupies hollow sites and does 
not change the symmetry of the outermost metal layers 
(Fig. [5|) . This situation results in strong interactions (see 
table ITT)) and corresponds to the incorporation of isolated 
or one-fold coordinated C atoms in the hollow sites, ob- 
served during our calculations. 

As the coverage exceeds 0.33 ML, the surface under- 
goes a reconstruction known as the "clock reconstruc- 
tion" described previously (Fig.[T]). This reconstruction. 
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TABLE 11: Binding energies (Ej) of a C atom on a Ni(lOO) slab. R is the distance between C and Ni atom. 



Efc (eV) R (A) 

Theory " -7.43, -6.50, -8.10, 1.91, 1.85, 1.77 

-7.69, -8.43 

Experiment " -7.35, -7.37, -7.55 1.80, 1.82, 1.89 

Present work -8.21 L90 

"From Ref. |3ll and references cited therein. 



which has also been studied theoretically [1^, [131 , is in- 
deed observed during our MC calculations (Fig.[3]). 

The adhesive properties of carbon clusters on Ni(lOO) 
surface then result from a competition between the for- 
mation of Ni-C and C-C bonds. For structures containing 
few C-C bonds, the competition between a carbide and a 
demixed system favors the formation of Ni-C bonds and 
the C atoms occupy the hollow sites of Ni(lOO) present 
at the surface. For planar flakes, C-C interactions pre- 
dominate and, as described above, the weak energy of 
adhesion is principally related to the energy gain of the 
(reconstructed) metallic surface below the flakes. Tube 
embryos have an intermediate behaviour: C-C interac- 
tions predominate on the caps and Ni-C anchoring bonds 
predominate at their edges. 

It is interesting to compare briefly our results with 
those obtained for the (111) surface. At low carbon con- 
centration the situation is similar, with predominant Ni- 
C interactions, even if, due to the different topology of 
the surfaces, the preferred position for the (111) surface 
for carbon is a subsurface octahedral site. At higher con- 
centration we expect the formation of a graphene sheet 
in quasi perfect epitaxial relationship with the Ni sub- 
strate and a vanishingly small energy of adhesion [2l| . 
We have also calculated the energy of adhesion of a small 
flat cluster made of 24 atoms. As in the case of the (100) 
surface we find a weak energy of adhesion equal to -0.26 
eV/C atom, and the local energy analyses show that the 



energy gain is concentrated on the edge atoms, without 
any important saturation of the dangling bonds. In their 
analysis of different caps on Ni (111) Reich et al. have 
also noticed that the relative energies involved are also 
fairly weak p^ . 



V. CONCLUSION 

In summary, we have presented tight-binding Monte 
Carlo simulations of the interaction of carbon structures 
with Ni(lOO) surfaces Our calculations show that during 
the first stages of the nucleation the incorporation of 
defects and the anchoring of carbon atoms to the surface 
are energetically favorable. Even in the presence of a 
metal substrate, the carbon dangling bonds parallel to 
the surface still exist and the metal-carbon system gains 
energy by reconstructing the underlying metallic surface, 
close to the well-known clock reconstruction. For curved 
structures, the formation Ni-C bonds is more efficient, 
and C atoms prefer to be located in semi-octahedral 
sites of Ni present at the surface. 

Part of this work has been supported by the Belgian 
Program on Interuniversity Attraction Poles (PAI5/1/1) 
on "Quantum Size Effects in Nanostructured Materials" . 



[1] C. Journet, W.K. Maser, P. Bernier, A. Loiseau, M. 
Lamy de la Chapelle, S. Lefrant, P. Deniard, R. Lee, 
and J. E. Fisher, Nature 388, 756 (1997). 

[2] A. Thess, R. Lee, P. Nikolaev, H. Dai, P. Petit, J. Robert, 
C. Xu, Y. H. Lee, S. G. Kim, A. G. Rinzler, D. T. Col- 
bert, G. E. Scuseria, D. Tomanek, J. Fischer, and R. E. 
Smalley, Science 273, 483 (1996). 

[3] J.- F. Colomer, G. Bister, L Willems, Z. Konya, A. Fon- 
seca, G. Van Tendeloo, and J. B. Nagy, Chem. Comm. 
14, 1343, (1999). 

[4] J. Gavillet, A. Loiseau, C. Journet, F. Willaime, F. 
Ducastelle, and J.-C. Charlier, Phys. Rev. Lett. 87, 
275504, (2001). 

[5] J. Gavillet, J. Thibault, O. Stephan, H. Amara, A. 
Loiseau, C. Bichara, J.-P. Gaspard and F. Ducastelle, 
J. Nanosci. Nanotech. 4, 346, (2004). 



[6] A. Maiti, C. J. Brabec, C. Roland, J. Bernholc, Phys. 

Rev. B 52, 14850, (1995). 
[7] A. Maiti, C. J. Brabec, J. Bernholc, Phys. Rev. B 55, 

6097, (1997). 

[8] Y. Shibuta, and S. Maruyama, Chem. Phys. Lett. 382, 
381, (2003). 

[9] F. Ding, A. Rosen, and K. Bolton, and , J. Chem. Phys. 

121, 2775, (2004); F. Ding, K. Bohon, and A. Rosen, J. 

Phys. Chem. B 108, 17369 (2004). 
[10] J. H. Onuferko, D. P. Woodruff, and B. W. Holland, Surf. 

Sci. 87, 357, (1979). 
[11] M. Bader, C. Ocal, B. Hillert, J. Haase, and A. M. Brad- 

shaw, Phys. Rev. B 35, 5900, (1987). 
[12] C. Klink, L. Olesen, F. Besenbacher, I. Stensgaard, E. 

Laegsgaard and N. D. Lang, Phys. Rev. Lett. 71, 4350, 

(1993). 



7 



[13] C. Klink, I. Stonsgaard, F. Besenbacher, and E. Laegs- 

gaard, Surf. Sci. 342, 250, (1995). 
[14] X. Fan, R. Buczko, A. A. Purctzky, D. B. Gcohcgan, ,]. 

Y. Howe, S. T. Pantelides, and S. J. Pennycook, Phys. 

Rev. Lett. 90, 145501, (2003). 
[15] J. Y. Raty, F. Gygi,, and G. Galli, Phys. Rev. Lett. 95, 

096103 (2005). 

[16] S. Reich, L. Li, and J. Robertson, Chem. Phys. Letters 

421, 469 (2006); S. Reich, L. Li, and J. Robertson, Phys. 

Rev. B 72, 165423, (2005). 
[17] V. L. Kuznetsov, A. N. Usoltseva, A. L. ChuviUn, E. D. 

Obraztsova, and J.-M. Bonard, Phys. Rev. B 64, 235401 

(2001). 

[18] H. Kanzow, and A. Ding, Phys. Rev. B 60, 11180 (1999); 

H. Kanzow, C. Lcnski, and A. Ding, Phys. Rev. B 63 

125402 (2001); J. Zhao, A Martinez-Limia, and P. B. 

Balbuena, Nanotechnology 16, S575 (2005). 
[19] Unfortunately it is difficult to reproduce the analysis 

from the data presented in the article. 
[20] H. Amara, These de Doctoral, Univ. Paris VI (2005); 

available at http://lem.onera.fr/Theses.html; H. 

Amara, ,J.- M. Roussol, C. Bichara, J.- P. Gaspard, and 

F. Ducastelle, to be published. 
[21] H. Amara, C. Bichara, and F. Ducastelle, Phys. Rev. B 



73, 113404, (2006). 
[22] N. Mctropohs, A. W. Rosenbluth, M. N. Roseubluth, 

A. Teller, and E. Teller, J. Chem. Phys. 21, 1087, (1953). 
[23] R. Haydock, V. Heine, and M.J. Kelly, J. Phys. C, 5, 

2845, (1972). 

[24] J. C. Shelton, H. R. Patil and J. M. Blakely, Surf. Sci. 

43, 493-520 (1974). 
[25] D. J. Siegel and J. C. Hamilton, Phys. Rev. B 68, 094105, 

(2003). 

[26] D. Alfc, S. de Gironcoh, and S. Baroni, Surf. Sci. 437, 
18, (1999). 

[27] S. Stolbov, S. Hong, A. Kara, and T. S. Rahman, Phys. 

Rev. B 72, 155423 (2005). 
[28] J. Fujii, G. Panaccione, I. Vobornik, G. Rossi, G. Tri- 

marchi and N. Binggch, Surf. Sci. 600, 3884, (2006). 
[29] H. S. Bengaard, J. K. N0rskov, J. Sehested, B. S. 

Clausen, L. P. Nielsen, A. M. Molenbroek, and J. R. 

Rostrup-Nielsen, J. Catal. 209, 365, (2002). 
[30] Y. H. Lee, S. G. Kim, and D. Tomanek, Phys. Rev. Lett. 

78, 2393, (1997). 
[31] Q. -M. Zhang, J. C. Wells, X. G. Gong, and Z. Zhang, 

Phys. Rev. B 69, 205413 (2004). 



